arXiv:1505.05803v3 [physics.bio-ph] 15 Nov 2015 


Particle diffusion in active fluids is non-monotonic in size 
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We experimentally investigate the effect of particle size on the motion of passive polystyrene spheres in suspensions of Es¬ 
cherichia coli. Using particles covering a range of sizes from 0.6 to 39 microns, we probe particle dynamics at both short and 
long time scales. In all cases, the particles exhibit super-diffusive ballistic behavior at short times before eventually transitioning 
to diffusive behavior. Surprisingly, we find a regime in which larger particles can diffuse faster than smaller particles: the particle 
long-time effective diffusivity exhibits a peak in particle size, which is a deviation from classical thermal diffusion. We also find 
that the active contribution to particle diffusion is controlled by a dimensionless parameter, the Peclet number. A minimal model 
qualitatively explains the existence of the effective diffusivity peak and its dependence on bacterial concentration. Our results 
have broad implications on characterizing active fiuids using concepts drawn from classical thermodynamics. 


1 Introduction 

The diffusion of molecules and particles in a fiuid is a process 
that permeates many^aspects of our lives including fog forma¬ 
tion in rain or snow®, cellular respiration®, and chemical dis¬ 
tillation processes®. At equilibrium, the diffusion of colloidal 
particles in a fiuid is driven by thermal motion and damped 
by viscous resistance®. In non-equilibrium systems, fiuctua- 
tions are no longer only thermal and the link between these 
fiuctuations and particle dynamics remain elusive®. Much ef¬ 
fort has been devoted to understanding particle dynamics in 
non-equilibrium systems, such as glassy materials and sheared 
granular mattei®. A non-equilibrium system of emerging in¬ 
terest is active matter. Active matter includes active fiuids, that 
is, fiuids that contain self-propelling parti cles, such as motile 
microorganisms®! catalytic colloids®ni and molecular mo- 
torsM These particles inject energy, generate mechanical 
stresses, and create flows w ithin the fiuid medium even in the 
absence of external forcingEGU^ Consequently, active fiuids 
display fascinating phenomena not seen in passive fiu ids, such 
as spontaneous flows® , ano malous shear viscosities®^, un¬ 
usual polymer swellingQ^G^, and enhanced fiuid mixingnSESl 
Active fiuids also play important roles in varied biological and 
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Fig. 1 (Color online) Trajectories of 2 particles (a) without 
bacteria and (b) with bacteria (c = 3 x 10^ cells/mL) for time 
interval 8 s. Trajectories of (c) 0.6 and (d) 16 ^um particles 
(c = 3 X 10^ cells/mL). Scale bar is 20 pm. 


























ecological settings, which include the contributions of suspen¬ 
sions of microorganisms to biofilm infectionsl^EII, biofouling 
of water-treatment systems^S, and biodegradation of environ¬ 
mental pollutants®!. 

The motion of passive particles in active fluids (e.g. suspen¬ 
sion of swimming microorganisms) can be used to investigate 
the non-equilibrium properties of such fiuids. For example, at 
short times particle displacement distributions can exhibit ex¬ 
tended non- Gaussi an tails as seen in experiments with algae 
C. reinhardti^^. At long times, particles exhibit enhanced 
diffusivities D^greater than Z)o, their thermal (Brownian) 
diffusivity These traits are a signature of the non¬ 

equilibrium nature of active fiuids; the deviation from equilib¬ 
rium also manifests in violations of the fiuctuation dissipation 
theorem®!. 

In bacterial suspensions, the enhanced diffusivity de¬ 
pends on the concentration c of bacteria. In their seminal 
work, Wu and Libchabei®! experimentally found that Dgff in- 
crease d linearly w ith c in suspensions of E. coli. Subsequent 
studies®®®^!! have observed that this scaling holds at low 
concentrations and in the absence of collective motion. In this 
regime, Dgff ca n be dec omposed into additive components as 
= Dq+ DAl 26 | 27 | 29| tll] are the thermal 

and active diffusivities, respectively. It has been proposed that 
the active diffusivity Da is a consequence of advection due to 
far-field i nterac tions with bacteria®! and may even be higher 

near walls®®! 

While a majority of studies have focused on the role of bac¬ 
terial concentration c on particle diffusion, the role of particle 
diameter d remains unclear. In the absence of bacteria, the 
diffusivity of a sphere follows the Stokes-Einstein relation. 
Do = ksT/fo, where ks is the Boltzmann constant, T is the 
temperature, and fo = 371 fid is the Stokes friction factoiE!in 
a fiuid of viscosity /i. In a bacterial suspension, this relation 
is no longer expected to be valid. Surprisingly, for large par¬ 
ticles (4.5 and 10 /im), Wu and Libchabei®! suggested that 
Deff scales as 1/J, as in passive fiuids. Recent theory and sim¬ 
ulation by Kasyap et al^ however do not support the l/d 
scaling and instead predict a non-trivial dependence of Dgff on 
particle size, including a peak in Dgff. This non-monotonic 
dependence of Dgff on particle size implies th at me asures of 
effective dif fusivi ties®! effective temperatures®®!, and mo¬ 
mentum fiux®®! intimately depend on the probe size and thus 
are not universal measures of activity. This has important im¬ 
plications for the common use of colloidal probes in gauging 
and characterizing t he act ivity of h ying m aterials, such as sus¬ 
pensions of bacteria®®!, biofilms®®!, and the cytoskeletal 
network inside cells®!, as well as in understanding transport 
in these biophysical setting. Despite the ubiquity of passive 
particles in active environments, the effects of size on particle 
dynamics in active fiuids has yet to be systematically investi¬ 
gated in experiments. 
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Fig. 2 (Color online) (a) Mean-square displacements (MSD) of 2 
jum particles versus time At for varying bacterial concentration c. 
Dashed line is a fit to Langevin dynamics (eqn (1)). (b) MSD for 
varying particle diameter d versus time for bacterial concentration 
c = 3.0 X 10^ cells/mL. The MSD peaks atd = 2 jum. The 
cross-over time T (arrows) increases with d. Solid lines are 
Langevin dynamics fits. 

In this manuscript, we experimentally investigate the effects 
of particle size d on the dynamics of passive particles in sus¬ 
pensions of Escherichia coli. Escherichia col^ are model 
organisms for bacterial studies and are rod-shaped cells with 
3 to 4 fiagella that bundle together as the cell swims forward 
at speed U approximately 10 /i/s. We change the particle size 
d from 0.6 /im to 39 /im, above and below the effective total 
length (L ^ 7.6 /im) of the E. coli body and flagellar bun¬ 
dle. We find that Dgff is non-monotonic in J, with a peak at 
2 < J < 10 /im; this non-monotonicity is unlike the previ¬ 
ously found \/d scaling®! and suggests that larger particles 
can diffuse faster than smaller particles in active fiuids. Fur¬ 
thermore, the existence and position of the peak can be tuned 
by varying the bacterial concentration c. The active diffusion 
Da = Deff — Do is also a non-monotonic function of d and can 
be collapsed into a master curve when rescaled by the quan¬ 
tity cUL!^ and plotted as a function of the Peclet number Pe 
= UL/D{) (cf Fig. 5(b)). This result suggests that the active 
contribution to particle diffusion can be encapsulated by an 
universal dimensionless dispersivity Da that is set by the ratio 





















of times for the particle to thermally diffuse a distance L and 
a bacterium to swim a distance L. 

2 Experimental Methods 

Active fluids are prepared by suspending spherical polystyrene 
particles and swimming E. coli (wild-type K12 MG1655) in 
a buffer solution (67 mM NaCl in water). The E. coli are 
prepared by growing the cells to saturation (10^ cells/mL) 
in culture media (LB broth, Sigma-Aldrich). The saturated 
culture is gently cleaned by centrifugation and resuspended 
in the buffer. The polystyrene particles (density p = 1.05 
g/cm^) are cleaned by centrifugation and then suspended in 
the buffer-bacterial suspension, with a small amount of sur¬ 
factant (Tween 20, 0.03% by volume). The particle volume 
fractions (j) are below 0.1% and thus considered dilute. The E. 
coli concentration c ranges from 0.75 to 7.5 xlO^ cells/mL. 
These concentrations are also considered dilute, correspond¬ 
ing to volume fractions 0 = cvb <1%, where Vb is the vol- 
umeE^of an E. coli cell body (1.4 jivo?). 

A 2 pi drop of the bacteria-particle suspen sion is stretched 
into a fluid film using an adjustable wire frameEGHEHlto a mea¬ 
sured thickness of approximately 100 pm. The film interfaces 
are stress-free, which minimizes velocity gradients transverse 
to the film. We do not observe any large scale collective be¬ 
havior in these Aims; the E. coli concentrations we use are 
below the values for which collective motion is typically ob- 
servedlSI 10^^ cells/ml ). Particles of different diameters 
(0.6 pm <d <39 pm) are imaged in a quasi two-dimensional 
slice (10 pm depth of focus). We consider the effects of parti¬ 
cle sedimentation, interface deformation, and confinement on 
particle diffusion and And that they do not significantly af¬ 
fect our measurements of effective diffusivity in the presence 
of bacteria (for more details, please see SI-§I). Images are 
taken at 30 frames per second using a lOX objective (NA 0.45) 
and a CCD camera (Sony XCDSX90), which is high enough 
to observe correlated motion of the particles in the presence 
of bacteria (Fig. 2) but small enough to resolve spatial dis¬ 
placements. Particles less than 2 pm in diameter are imaged 
with fluorescence microscopy (red, 589 nm) to clearly visu¬ 
alize particles distinct from E. coli (2 pm long). We obtain 
the particle positions in tw o dimensions over time using par¬ 
ticle tracking methods^ESI All experiments are performed at 
To = 22°C. 

3 Results and Discussion 

3.1 Mean Square Displacements 

Representative trajectories of passive particles in the absence 
and presence of E. coli are shown in Fig. 1 for a time inter¬ 
val of 8 s. By comparing Fig. 1(a) (no E. coli) to Fig. 1(b) 


(c = 3 X 10^ cells/mL), we readily observe that the presence 
of bacteria enhances the magnitude of particle displacements 
compared to thermal equilibrium. Next, we compare sample 
trajectories of passive particles of different sizes d, below and 
above the E. coli total length L ^ 7.6 pm. Figure 1(c) and 1(d) 
show the magnitude of particle displacement for J = 0.6 pm 
and J = 16 pm, respectively. Surprisingly, we And that the 
particle mean square displacements in the E. coli suspension 
are relatively similar for the two particle sizes even though the 
thermal diffusivity Dq of the 0.6 pm particle is 35 times larger 
than that of the 16 pm particle. The 16 pm particles also 
appear to be correlated for longer times than the 0.6 pm parti¬ 
cles. These observations point to a non-trivial dependence of 
particle diffusivity on d and are illustrated in sample videos 
included in the Supplementary Materials. 

To quantify the above observations, we measure the mean- 
squared displacement (MSD) of the passive particles for vary¬ 
ing E. coli concentration c (Fig. 2(a)) and particle size d (Fig. 
2(b)). Here, we deflne the mean-squared particle displacement 
as MSD(Ar) = (|r(rR +Ar) — r(rR)p), where the brackets de¬ 
note an ensemble average over particles and reference times 
^R. For a particle executing a random walk in two dimensions, 
the MSD exhibits a characteristic cross-over time T, corre¬ 
sponding to the transition from an initially ballistic regime for 
Ar <C T to a diffusive regime with MSD ~ 4Z)eff Ar for ArT. 

Figure 2(a) shows the MSD data for 2 pm particles at vary¬ 
ing bacterial concentrations c. In the absence of bacteria (c = 0 
cells/mL), the fluid is at equilibrium and Dgff = Dq. For this 
case, we are unable to capture the crossover from ballistic 
to diffusive dynamics due to the lack of resolution: for col¬ 
loidal particles in water, for example, cross-over times are 
on the order of nanoseconds and challenging to measureESI. 
Experimentally, the dynamics of passive particles at equilib¬ 
rium are thus generally diffusive at all observable time scales. 
We lit the MSD data for the J = 2 pm case (with no bacte¬ 
ria) to the expression MSD = 4 Z)oAl and And that Dq 0.2 
pm^/s. This matches the theoretically predicted value from 
the Stokes-Einstein relation; the agreement (V) can be visu¬ 
ally inspected in Fig. 3(a). 

In the presence of E. coli, the MSD curves exhibit a ballis¬ 
tic to diffusive transition, and we And that the cross-over time 
T increases with c. For Ar T, the MSD ^ Ar with a long¬ 
time slope that increases with bacteria concentration c. Ad¬ 
ditionally, the distribution of particle displacements follows a 
Gaussian distribution (see SI-§II A for details and measures of 
the non-Gaussian parameter). These features, MSD ^ Ar and 
Gaussian displacements, indicate that the long-time dynam¬ 
ics of the particles in the presence of E. coli is diffusive and 
can be captured by a physically meaningful effective diffusion 
coefficient Dgff. 

We next turn our attention to the effects of particle size. For 
varying particle diameter J at a flxed bacterial concentration 
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Fig. 3 (Color online) (a) Effective particle diffusivities Dgff versus particle diameter d at varying c. The dashed line is particle thermal 
diffusivity Dq. (b) The crossover-time T increases with d, scaling as approximately J", where \/2<n< 1. This is not the scaling in passive 
fluidsl^ where T oc d^. 


(c = 3 X 10^ cells/ml), the MSD curves also exhibit a ballistic 
to diffusive transition, as shown in Fig. 2(b). Surprisingly, we 
find a non-monotonic behavior with d. For example, the MSD 
curve for the 2 /im case sits higher than the 39 /im case and the 
0.6 /im case. This trend is not consistent with classical diffu¬ 
sion in which MSD curves are expected to decrease monoton- 
ically with d (Dq \/d). We also observe that the cross-over 
time T increases monotonically with d. As we will discuss 
later in the manuscript, the cross-over time scaling with d also 
deviates from classical diffusion. 

3.2 Diffusivity and Cross-over Times 

We now estimate the effective diffusivities Dgff and cross-over 
times T of the passive particles in our bacterial suspensions. 
To obtain Dgff and T, we fit the MSD data shown in Fig. 2 to 
the MSD expression attained from the generalized Langevin 
equation® that is 

MSD(Ar)=4£)effAr(l-^(l-e“^)). (1) 

Equation 1 has been used previously to interpret the diffu¬ 
sion of bacteria® as well as the diffusion of particles in 
films with bacteria®. In the limit of zero bacterial con¬ 
centration, Da = 0 and eqn (1) reduces to the formal solu¬ 
tion to the Langevin equation for passive fiuids, MSD(Ar) = 
4DoAr (l — ^ (l — with To = t(c = 0). For more 

details on the choice of our model, see SI-§III. 

Figure 3 shows the long-time particle diffusivity Dgff (Fig. 
3(a)) and the cross-over time T (Fig. 3(b)) as a function of d 
for bacterial concentrations c = 0.75, 1.5, 3.0 and 7.5 x 10^ 
cells/mL. We find that for all values of d and c considered 
here, Dgff is larger than the Stokes-Einstein values Dq at equi¬ 
librium (dashed line). For the smallest particle diameter case 


(J = 0.6 /im), Deff nearly matches Dq. This suggests that 
activity-enhanced transport of small (J < 0.6 fim) particles 
or molecules such as oxygen, a nutrient for E. coli, may be 
entirely negligible®. For more information, including figures 
illustrating the dependence of Dgff on c and comparisons be¬ 
tween our measured effective diffusivities and previous exper¬ 
imental work, see sections SI-§IIB and SI-§IIC in the supple¬ 
mental materials. 

Figure 3(a) also reveals a striking feature, a peak Dgff in d. 
Our data demonstrates that, remarkably, larger particles can 
diffuse faster than smaller particles in suspensions of bacteria. 
For example, at c = 7.5 x 10^ cells/mL (Q) the 2 fim parti¬ 
cle has an effective diffusivity of approximately 2.0 jdm^/s, 
which is nearly twice as high as the effective diffusivity of the 
0.86 /im particle, Dgff =1.3 /im^/s. We also observe that the 
peak vanishes as c decreases. In fact, for the lowest bacterial 
concentration (c = 0.75 x 10^ cells/mL), there is no peak: Dgff 
decreases monotonically with d. Clearly, Dgff does not scale 
as 1/d. 

Figure 3(b) shows the cross-over times T characterizing the 
transition from ballistic to diffusive regimes as a function of 
particle size d for varying c. We find that the values of T in¬ 
crease with d and c. We note that the variation of T with c (Si- 
Fig. 3(b)) does not follow a linear form. Instead, the data sug¬ 
gests possible saturation of T for suspensions of higher - but 
still dilute - concentrations. The cross-over time T scales with 
particle diameter approximately as T ^ d^, where 
This cross-over time does not correspond to the inertial relax¬ 
ation of the particle. Therefore, this scaling does not follow 
the trend seen for passive particles at thermal equilibrium®, 
where T, being the Stokes relaxation time (order 1 ns), scales 
as m//o oc with particle mass m oc 

highlights that in active fiuids the super-diffusive motion of 
the passive particles cannot be ignored - even for time scales 
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Fig. 4 (Color online) (a) Distribution of 2 ^um particle speeds p{v) follow a Maxwell-Boltzmann distribution (solid curves) with clear peaks 
that shift right as c increases, (b) The effective temperature Teff extracted by fitting p{v) data to eqn (2) (A) matches those obtained from an 
extended Stokes-Einstein relation (Q)- (c) The effective temperature Tgff increases with d for varying c. 


as large as a second - and that the time scales over which 
diffusive motion is valid (Ar > t) depends on the size d of 
the particle. Further implications of a particle size-dependent 
cross-over time will be discussed below. 

3.3 Effective Temperature 

Our data so far suggests that particle dynamics in bacterial sus¬ 
pensions, while having an anomalous size-dependence (Figs. 
1,2,3), maintain the characteristic super-diffusive to diffusive 
dynamics for passive fluids® The long-time diffusive behav¬ 
ior (Fig. 2) and enhancement in (Fig. 3(a)), which is 
rooted in particle-bacteria encounters, suggest that the parti¬ 
cles behave as if they are suspended in a fluid with an effec¬ 
tively higher temperature. 

In order to further explore the concept of effective temper¬ 
ature in bacterial suspensions, we measure the distribution of 
particle speeds p{v) as a function of bacterial concentration 
c. The particle speed distributions determine the mean kinetic 
energy of the particles. If the distribution follows a Maxwell 
Boltzmann form - as is always the case in fluids at equilibrium 
- the mean kinetic energy is related to the thermodynamic 
temperature via the equipartition theorem. Such a relationship 
may not always exist for out-of-equilibrium fluids. 

Figure 4(a) shows p{v) for J = 2 fim case for a range of c. 
We define the particle speeds v over a time interval of 0.5 s. 
This time interval is greater than the crossover times for the 
2 /im particles (cf. Si-Fig. 3) to ensure that a particle sam¬ 
ples multiple interactions with bacteria and exhibits diffusive 
behavior. In the absence of bacteria, the system is in thermal 
equilibrium and p{v) follows the two dimensional Maxwell- 
Boltzmann distribution, 

2 

—mv^ 

p{y) = vm(^5reff)“^e^^57eff , (2) 


with peak speeds Vmax = where m is the mass of 

the polystyrene particle. Fitting the p{v) data in the absence 
of bacteria (V in Fig. 4(a)) to eqn (2) yields Tgff ^ Tq, as ex¬ 
pected. 

In the presence of bacteria, the particle speed distributions 
also follow a Maxwell-Boltzmann form, eqn (2), with peaks 
that shift toward higher values of v as the E. coli concentration 
c increases. Because our data follows the Maxwell-Boltzmann 
form (Fig. 4(a)) for all c, this indicates that there is no corre¬ 
lated motion at long times as is the case in swarming bacterial 
suspensions, for which the particle speed distribtuions exhibits 
an exponential decay®. We also note that the power spectra of 
particle speeds (SI-§II D) are reasonably flat, consistent with 
white-noise forcing and an absence of correlated motion. Fig¬ 
ure 4(a) thus indicates that an effective temperature Tgff can be 
defined from p{y) and is increasing with the bacterial concen¬ 
tration. 

To quantify this ‘enhanced’ temperature and the deviation 
from equilibrium behavior, we fit the p{v) data (Fig. 4(a)) to 
eqn (2). These fits allow us to obtain Tgff as a function of c, 
as shown in Fig. 4(b) (A). Also shown in Fig. 4(b) are the 
values extracted of Teff (Q) an extended Stokes-Einstein 
relation, = ?>niidDQfflkB, where Z)eff are the values from 
the MSD fits shown in Fig. 3(a) and jl is the viscosity of the 
solvent (/i ~ 1 mPa s). We find that both estimates of T^ff are 
higher than room temperature Tq and increase linearly with c. 

For a fluid at equilibrium, the temperatures from these two 
methods, MSD and p{v), are expected t o be the same but not 
for systems that are out of equilibriumSEIElES, such as the 
bacterial suspensions investigated here. The good agreement 
in Fig. 4(b) suggests that Tgff may be a useful signature of 
bacterial activity with some analogies to equilibrium systems. 
Note that in defining the effective temperature via the general¬ 
ized Stokes-Einstein relation, we have assumed the viscosity 
to be a constant and independent of bacterial concentration. 
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which may not be true when the bacterial concentrations are {E. coli) as well as challenges in gauging activity using passive 
sufficiently high®!. Figure 4(b) suggests that an unchanging particles, 
viscosity is a valid assumption for our system. 
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Fig. 5 (Color online) (a) Active diffusivities — Dq are 

non-monotonic with particle size for varying concentrations of 
bacteria, (b) Scaled hydrodynamic diffusivity Da = Dp^/cUL^ 
collapses with Peclet number Pe = C/L/Dq. The maximum Da 
occurs at Pca between 450 and 4500. For 100 < Pe < Pca, Da 
scales as Pe^, where a'^1. 


Next, we investigate the role of particle diameter in the 
suspension effective temperature. Figure 4(c) show the val¬ 
ues of Deff estimated from an extended Stokes-Einstein re¬ 
lation as a function of d for different values of c. Surpris¬ 
ingly, we find that Dgff increases with particle size J, which is 
clearly different from thermally equilibrated systems where 
temperature does not depend on the probe size. We note 
that for the largest particle diameter, J = 39 fim, estimated 
Deff values are approximately 100 times greater than room 
temperature Tq = 295 K, consistent with previous reports!^. 
The dependence of Dgff on particle size d (Fig. 5(a)) may 
be understood through the extended Stokes-Einstein relation, 
which yields T^ff = D^ff/o/ks = Dq -f {?>7tilDp^/kB)d, where 
Tq = {37riiDo/kB)d. If Da were independent of J, then Dgff 
would be linear in d. However, a linear fit does not adequately 
capture the trend. This hints at a particle size dependent active 
diffusivity Da, which has been predicted in recent theoryEZI. 
This variation of Dgff with d highlights the interplay between 
particle size and the properties of the self-propelling particles 


3.4 Active Diffusivity of Passive Particles in Bacterial 
Suspensions 

To explore the aforementioned dependence of Da on particle 
size, we plot Da = Dgff — Dq with d. Here, Dq is the parti¬ 
cle thermal diffusivity (in the absence of bacteria) and is ob¬ 
tained from the Stokes-Einstein relation Dq = . Indeed, 

as shown in Eig. 5(a), Da exhibits a non-monotonic depen¬ 
dence on d for all c. 

To understand the observed dependence of the active diffu¬ 
sivity on the particle size, we consider the relevant time scales 
in our particle/bacteria suspensions, namely: (i) the time for 
the particle to thermally diffuse a distance L equal to the total 
bacterial length l} (ii) the time for a bacterium to swim at 
a speed D for a distance L given by L/D, and (iii) the mean run 
time which is the inverse of the tumbling frequency 0)^^. Di¬ 
mensionless analysis then suggests there are two independent 
time parameters: (i) the ratio of the first two above, which is 
the Peclet number, Pe = DL/Dq, and (ii) T* = D/ (OrL, which 
is the ratio of the run length D/(0r to the bacterial length. 

It is these two parameters, Pe and T*, that govern the parti¬ 
cle dynamics in our experiments. When the Peclet number is 
much less than one, then the thermal particle diffusion dom¬ 
inates and transport by the bacteria is ineffective. When the 
Peclet number is much larger than one, then thermal diffusion 
is negligible and the transport is due to the convection from 
bacteria.The swimming speed D, tumble frequencies , 
and combined length of the cell bod y and fiagellar bundle L 
are estimated from prior experiments^^^ with E. coli as ap¬ 
proximately 10 /im/s, 1 s“^ and L « 7.6 /im, respectively. 
Thus, in our experiments, the Peclet number varies from ap¬ 
proximately 130 to 8600, via the particle bare diffusivity Dq 
(through the particle diameter). We note that one stain of bac¬ 
teria is used - thus, the run length and bacteria size do not 
change in our experiments, and consequently, T* ^ 1.8 is a 
constant. 

In order to gain insight into the non-trivial dependence of 
Da on J, we scale out the concentration dependence by in¬ 
troducing the dimensionless active diffusivity Da = Da/cUL^ 
and plot it against the Peclet number, Pe. Eigure 5(b) shows 
that all the active diffusion Da versus d data shown in Eigure 
5(a) collapses into a single master curve, thereby indicating 
that Da is independent of c, at least for dilute suspensions in¬ 
vestigated here. 

We find that for Pe < 10^, the values of Da initially in¬ 
creases with increasing Pe (or particle size) and follows a scal¬ 
ing Da ^ Pe^. The observed increase of Da with Pe may be 
due to the decreasing particle Brownian motion, which allows 
the particle’s motion to be correlated with the bacterial veloc- 


















ity disturbances for longer times. 

In the limit of Pe ^ oo, the particle’s Brownian motion van¬ 
ishes and the particle displacements are dominated by the con¬ 
vective transport via bacteria-particle interactions. Thus Da is 
expected to be independent of Pe and depend only on the pa¬ 
rameter T*, defined here as the ratio of the run length U/ (Dr io 
the bacterial length L. Given that our experiments correspond 
to a constant T* ^ 1.8, our experiment ally measured Da agree 
well with recent theoretical predictions^ for very large Peclet 
(as Pe ^ oo) for t* = 1 and T* = 4, as shown in Fig. 5(b). An 
increase in the run time (or T*) would increase the asymptotic 
value of the scaled active diffusivity Da (see SI-§IV). 

An important feature of the data shown in Fig. 5(b) is the 
peak in Da at PeA ~ 10^. The existence of such a peak has 
been predicted in a recent theory/simulation investigational^ 
and our data agrees at least qualitatively with such predic¬ 
tions. We note that the predicted peak in Da happens at a 
Pe ~ 0(10), which is smaller than our experimental values of 
(Pe ^10^). The appearance of the peak in Da in our data may 
be due to the weak but non-zero effects of Brownian motion, 
which allows particles to sample the bacterial velocity field in 
such a way that the mean square particle displacements and 
correlation times are higher compared to both the very high 
Pe (negligible Brownian motion) as well as the very low Pe 
(Brownian dominated)!^!. This feature (i.e. peak in Da) is sur¬ 
prising because it suggests an optimum particle size for max¬ 
imum particle diffusivity that is coupled to the activity of the 
bacteria. 


D^fi = Da + {cL^){UL)DA. (3) 

The criterion for the existence of a maximum Dgff is obtained 
by taking the derivative of eqn (3) with respect to the Peclet 
number and setting the derivative to zero. In order to estimate 
Da (and its slope), we fit the data in Fig. 5(b) near the peak 
in the range 200 < Pe < 4000 with a second order polynomial 
equation. We find that a peak exists in Dgff if 


cL^ > 0.4. (4) 

For the bacterial length used here L = 7.6 /im, this yields c ^ 
0.9 X 10^ cells/mL, which is in quantitative agreement with 
the concentration range (0.75 x 10^ cells/mL < c < 1.5 x 10^ 
cells/mL) in which the peak in Dgff emerges in our data (Fig. 
3(a)). As described in SI-§V, we find that the location of the 
Z)eff P^ak in d here defined as is given by 


jmax , 
"eff ' 


!(/a 


1 


5cL3 - 2 


(5) 


where Ja corresponds to the Peclet number PeA ~ 1000 at 
which Da is maximum. For the E. coli used here, Ja = 
^rPeA/3;r/if/L 6 jim. Note that d^^^ is always less than 
or equal to Ja and increases with c, consistent with our exper¬ 
imental observations (Fig. 3(a)). In general the criterion, eqn 
(4), and d^f^, eqn (5), will depend on T*. For suspensions of 
bacteria, the universal curve of Da informs when and where a 
peak in the particle diffusivity occurs. 


4 Maximum Particle Effective Diffusivity D^ff 

Our data (Fig. 5(b)) shows that the dimensionless active dif¬ 
fusion Da collapses unto a universal curve with a peak in Pe 
for all bacterial concentrations. Our data also shows that, with 
the exception of the lowest bacterial concentration c, the parti¬ 
cle effective diffusivity D^ff exhibits a peak in J, which varies 
with c (Fig. 3(a)). For instance, for c = 1.5 x 10^ cells/mL, 
the peak is at approximately 2 /im, while for larger concentra¬ 
tions (c = 1.5 X 10^ cells/mL), the peak (as obtained by fitting 
the data to a continuous function) shifts to higher values of 
d. This suggests that one can select the particle size which 
diffuses the most by tuning the bacteria concentration. 

In what follows (see also SI-§V), we provide a prediction, 
based on our experimentally-measured universal curve of Da 
with Pe (Fig. 5(b)) for the existence as well as the location 
of the peak of D^ff in d. As noted before, the particle effec¬ 
tive diffusivity D^ff can be described as the linear sum of the 
particle thermal diffusivity Do, which is independent of c and 
decreases with d, and the active diffusivity Da, which is linear 
in c and non-monotonic in d, through the particle-size depen¬ 
dent Da. Therefore, we can recast the effective diffusivity as 


5 Conclusions 

In summary, we find that the effective particle diffusivity D^ff 
and temperature T^ff in suspensions of E. coli show strong de¬ 
viations from classical Brownian motion in the way they de¬ 
pend on particle size d. For example. Fig. 3(a) shows that 
Z)eff depends non-monotonically in d and includes a regime 
in which larger particles can diffuse faster than smaller parti¬ 
cles. The existence as well as the position of a Dgff P^ak in d 
can be tuned by varying the bacterial concentration c. We also 
find that the cross-over time T increases with particle size and 
scales as approximately d^, where l/2<^< 1, as shown in 
Fig. 3(b). 

Measures of Tgff obtained from either an extended Stokes- 
Einstein relation or particle speed distributions seem to agree 
quite well (Fig. 4(b)). This is surprising since this kind of 
agreement is only expected for systems at equilibrium. The 
good agreement between the measurements suggests that Tgff 
may be a useful signature of bacterial activity. However, un¬ 
like thermally equilibrated systems, Tgff varies with size d 
(Fig. 4(c)). This non-trivial dependence of both D^a and Tgff 
on particle size d implies that these common gauges of activ¬ 
ity are not universal measures. Nevertheless, our data suggest 







that one can define optimal colloidal probes of activity for sus¬ 
pensions of bacteria, which correspond to Pe = f/L/Z)o ~ 10^. 
At these Pe values, 5 a is maximized, which provides ample 
dynamical range (magnitude of the signal). Also at these Pe 
values, the cross-over time T is still relatively small, which al¬ 
lows for adequate temporal resolution. Both of these features 
are important in using passive particles to characterize a spa¬ 
tially and temporally varying level of activity in materials. 

Our anomalous particle-size dependent results in active fiu- 
ids has important implications for particle sorting in microfiu- 
idic devices, drug delivery to combat microbial infections, re¬ 
suspension of impurities and the carbon cycle in geophysical 
settings populated by microorganisms. A natural next step 
would be to study the role of external fields such as gravity 
or shear in influencing particle transport in these active envi¬ 
ronments. 
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